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ABSTRACT 

Context. The transient X-ray source XMMU J004215. 8+41 1924 within M31 was found to be in outburst again in the 2010 May 27 
Chandra observation. We present results from our four Chandra and seven Swift observations that covered this outburst. 
Aims. X-ray transient behaviour is generally caused by one of two things: mass accretion from a high mass companion during some 
restricted phase range in the orbital cycle, or disc instability in a low mass system. We aim to exploit Einstein, HST, Chandra and 
Swift observations to determine the nature of XMMU J004215.8+41 1924. 

Methods. We model the 2010 May spectrum, and use the results to convert from intensity to counts in the fainter Chandra observations, 
as well as the Swift observations; these data are used to create a lightcurve. We also estimate the flux in the 1979 January 13 Einstein 
observation. Additionally, we search for an optical counterpart in HST data. 

Results. Our best X-ray positions from the 2006 and 2010 outbursts are 0.3" apart, and 1.6" from the Einstein source; these outbursts 
are likely to come from the same star system. We see no evidence for an optical counterpart with itib S25.5; this new limit is 3.5 
magnitudes fainter than the existing one. Furthermore, we see no V band counterpart with m v 5=26. The local absorption is ~7 times 
higher than the Galactic line-of-sight, and provides ~2 magnitudes of extinction in the V band. Hence M\ S -0.5. Fits to the X-ray 
emission spectrum suggest a black hole primary. 

Conclusions. We find that XMMU J004215.8+411924 is most likely to be a transient LMXB, rather than a HMXB as originaly 
proposed. The nature of the primary is unclear, although we argue that a black hole is likely. 
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1. Introduction 

The bulge region of M31, the nearest spiral galaxy neighbour, is 
one of the best laboratories in the Universe for studying X-ray 
binaries. Accordingly, it has been observed hundreds of times 
by various X-ray observatories over the past 30 years. In the last 

. ~10 years alone, it has been observed 120 times by Chandra, 
90 times by Swift and 31 times by XMM-Newton. Most of the 
Chandra and XMM-Newton obbservations were short, monitor- 
ing observations looking for transient X-ray sources, while the 
Swift observations general ly followed these transients (see e.g. 

! IWilliams et alj|2rj05ll200rjh . 

XMMU J004215.8+411924 was identified as a new X-ray 
trans ient in the 2006, August 9 observation of M3 1 (|Haberl et all 
120061) . with a positional uncertainty of 2". lHaberl et al.l (12006) 
found that the emission spectrum was well described by a power 
law with photon index 1.57, suffering absorption equivalent to 
4.2xl0 21 H atom cm ; the unabsorbed 0.5-10 keV luminos- 
ity was 9.1xl0 37 erg s" 1 , assuming a distance of 780 kpc. The 
follow-up Swift observation made on 2006 September 1 re- 
vealed a UV count erpart within the X-ray error circle, leading 
lHaberl et alJ d2006h to identify XMMU J004215.8+41 1924 as a 
high mass X-ray binary (HMXB). A further Swift observation 
on 2006 September 1 1 yielded only thre e photons from XMMU 
J004215.8+41 1924 dPietsch et al.ll2006h : this corresponded to a 
0.5-10 keV luminosity of < 5 x 10 36 erg s _I when assuming the 
above emission model. 



iGalache et alJ (]2006) subsequently reported a Chandra de- 
tection of XMMU J004215.8+411924in a July 31 observation, 
finding the 0.9-6 keV spectrum to be well modeled by a power 
law model with photon index ~ 1 .8, with absorption equivalent to 
4.4xl0 21 H atom cnT 2 . They found the 0.5-10 keV luminosity 
to be l.lxlO 38 erg s" 1 . 

IVoss et all d2008l) examined the Chandra, Swift and XMM- 
Newton observations of three transients in M31, including 
XMMU J004215.8+411924. The June 2 XMM-Newton obser- 
vation made no firm detection of J004215.8+41 1924; hence 
IVoss et all (120081) constrainted the outburst duration to 40-79 
days. They refined the source positi on by registering the C handra 
data with the 2MASS catalogue of lSkrutskie et all d2006h . They 
obtained RA(J2000) = 00:42:16.1, Dec(J2000) = +41:19:26.7, 
with a lcr uncertainty of 0.5". This new position w as ~ 4" 
from the c o unte rpart identified bv lHaberl et alJ (|2006), allowing 
IVoss et all (120081) to reject this association. They searched for 
a counterpa rt in the local group galaxy survey (LGGS) images 
provided by iMassev et al.l (2006), and found nothing with V 5= 
22; using a distance modulus of 24.46 and 0.4 mag of extinction, 
they could not rule out a Be companion star. 

In iNooraee et al.l d2010h we reported a new outburst in the 
2010 May 27 Chandra observation within 0.5" of the position 
of XMMU J004215.8+411924, and identified it as a recurring 
ttansient. In this paper we present detailed analysis of the seven 
Swift and four Chandra observations made between 2010 March 
5 and 2010 July 20. We also search for a counterpart in the 2006 
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Table 1. Journal of observations. For each observation we give 
the instrument, observation number, date, exposure time and 
number of net source photons 



Instrument 


Obs 


Date 


Exposure 


Net counts 


Einstein HRI 


579 


1979-01-13 


29 ks 


29 


Chandra ACIS-I 


7139 


2006-07-31 


4ks 


433 


Chandra ACIS-I 


11279 


2010-03-05 


4ks 


5 


Chandra ACIS-I 


11838 


2010-05-27 


4ks 


279 


Swift XRT 


0031255012 


2010-06-06 


4ks 


44 


Swift XRT 


0031255013 


2010-06-09 


4ks 


37 


Swift XRT 


0031255014 


2010-06-12 


4ks 


24 


Swift XRT 


0031255015 


2010-06-15 


1.7 ks 


9 


Swift XRT 


0031255016 


2010-06-18 


3ks 


6" 


Chandra ACIS-I 


11839 


2010-06-23 


4ks 


27 


Swift XRT 


0031255018 


2010-06-24 


4ks 


1" 


Chandra ACIS-I 


11840 


2010-07-20 


4ks 


~0 fl 



8 Not a secure detection. Derived luminosities are 3cr upper limits. 



August 27 HST observation, and present evidence for a previous 
outburst in 1979 detected with Einstein. 

We discuss the observations and data analysis in Section 2. 
We then present our results in Section 3 and discuss in Section 
4 whether XMMU J004215. 8+41 1924 is a HMXB displaying 
orbital variability, or a transient low mass X-ray binary (LMXB) 
with unstable disc accretion. 

2. Observations and data analysis 

A journal of observations is provided in Table Q] 

2.1. Analysis of Chandra data 

The 2006 July 31 Chandra observation (7139) o f XMMU 
J004215. 8+41 192 has already been discussed by IVoss et all 
(2008). However, we registered the image to the LGGS B band 
image of M31 Field 6, which has positional uncertainties of 
0.25", using g lobular clusters (GCs) in the Revise d Bologna 
Catalogue V4 dGalleti et alJl2004l l2006l 120071 120091) that were 
X-ray bright. To do this, we used the iraf tool imcentroid to de- 
termine the position of each cluster in the Chandra and LGGS 
observations in image coordinates, then checked the positional 
uncertainties of each GC. We then used xy2sky to get the sky 
coordinates for each GC in the X-ray and LGGS images. This 
allowed us to map the LGGS coordinates to the X-ray positions 
with ccmap. The position of XMMU J004215.8+41 192 was de- 
termined in the corrected image using imcentroid, and the posi- 
tional uncertainties were combined with the 0.25" of the LGGS 
coordinate system. 

For the 2010 May 27 Chandra observation (11838), we lo- 
cated XMMU J004215. 8+41 192 as above. In addition, we ex- 
tracted source and background spectra, with related products, 
using ciao ver 4.2. These spectra were analysed with xspec ver 
12.6. XMMU J004215.8+41 192 was too faint for spectral mod- 
elling in the 2010 March 5 , 2010 June 24 and 2010, July 20 
observations; instead, we derived a conversion from intensity to 
flux for each observation using the best fit emission model from 
Obs 1 1838, as described below. 

2.2. Analysis of Swift data 

We were awarded seven observations of M3 1 with Swift dur- 
ing 2010, June at three day intervals to monito r the transient 
CXOM31 J004253. 1+4122 dHenze et all 12009): results from 



that transient will be presented by Nooraee et al. (in prep). This 
work uses data from the X-ray Telescope (XRT), taken in pho- 
ton counting (pc) mode. Unfortunately, the 2010 June 21 ob- 
servation had an exposure time of just 35s in pc mode, so this 
observation was ignored. XMMU J004215.8+41 1924 was at a 
large off-axis angle that varied between observations (~7— 10'). 

The Swift data were analysed with xselect version 2.4a. We 
were unable to use the recommended source extraction radius 
of 47" (containing ~ 90% of the source photons) due to crowd- 
ing. Instead we used a circular source extraction region with ra- 
dius 20", and a concentric, annular background region with in- 
ner radius 20" and outer radius 28.28". Lightcurve analysis was 
not productive, since the Swift data are obtained in chunks of a 
few hundred seconds, and very few photons were collected each 
time. 

Since there were too few photons for spectral analysis 
in each Swift observation, we converted from background- 
subtracted intensity to 0.5-10 keV flux by calculating the flux 
equivalent to 1 count s _1 assuming the best fit model derived 
from our analysis of the Obs 1 1838 Chandra data. To do this we 
obtained source and background spectra for each observation, 
and created an ancillary response file using the tool xrtmkarf; 
the appropriate canned response matrix was used. For each ob- 
servation we loaded these files into xspec version 12.6 and nor- 
malised the best fit emission model to give 1 count s _1 . This 
enabled us to obtain the unabsorbed 0.5-10 keV flux equivalent 
to 1 count s _1 , which we refer to as the conversion factor. 

2.3. Analysis of HST data 

We triggered two HST observations after the 2006 out- 
burst, designed to observed the counterpart in the on phase 
phase and off phase. The first observation was made with 
the ACIS/WFC on 2006 August 27 using the F435W filter. 
Unfortunately, the ACS was not operational during our second 
observation in 2007 July, and the observation was made with 
WFPC/FIX-1 using the F439W filter. However, the position of 
XMMU J004215.8+411924 was serendipitously covered by an 
ACS/WFC1 observation in 2004, October using the F555W fil- 
ter. 

We obtained the drizzled images of these observations from 
the Hubble legacy archive, and registered them with the B image 
of M31 Field 6 in the LGGS; we used stars that were visible in 
both images, but not so bright that the centroid determination 
had large uncertainties. 

3. Results 

3.1. Fitting the 2010 May Chandra spectrum 

The 2010 May observation of XMMU J004215.8+41 1924 
yielded 279 net source counts. We grouped the source spectrum 
to get a minimum of 15 counts per bin. The best fit power law 
emission model had a photon index of 1.8+0.5, with absorption 
equivalent to 6+3 x 10 21 H atom crrT 2 ;;^ 2 /dof =7/15. These pa- 
ramete rs are entirely consistent with those found by IVoss et all 
(120081) . 

Th e XMM-Newton observation analysed by IVoss et all 
(2008) yielded the tightest constraints on the absorption 
(4.2+0.5 x 10 21 H atom cm -2 ); we therefore fixed the absorp- 
tion to 4.2xl0 21 H atom cirT 2 for our Chandra spectrum. This 
gave a best fit photon index of 1.6+0.2, and^ 2 /dof =8/16. This 
yielded an unabsorbed 0.5-10 keV luminosity of 1.2+0.3 x 10 38 
erg s _1 . This fit is presented in Fig.Q] 
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Fig. 1. Best fit power law model to the Obs 1 1838 Chandra spec- 
trum, with residuals; ^- 2 /dof = 7/16. A^h was fixed to 4.2xl0 21 
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Fig. 2. 0.3-10 keV lightcurve of XMMU J004215.8+41 1924 
from Swift (filled circles) and Chandra (hollow circles). Arrows 
represent 3<x upper limits. The y axis is log scaled for clarity. 



We also fitted the spectrum with a disk blackbody model, 
since black holes in outburst are ofte n in the high soft state, 
with t hermal emission spectra (see e.g. McClinto ck & Remillardl 
2006, and references within). The best fit spectrum had an inner 
disc temperature of 1.4+0.3 keV; Nu was fixed to4.2xl0 21 atom 
cirT 2 as before; ^f 2 /dof = 10/16. This temperature is somewhat 
high f or a high state black hole binary (McClinto ck & Rem illard 
2006), and the fit is rather worse than the power law fit. Hence 
we prefer the power law fit, but cannot rule out a thermal fit. 

A neutron star accreting at ~ 10 38 erg s _1 has a two com- 
ponent emission spectrum, often characterised as a blackbody 
and a power law, with t he blackbody contributing ~ 10-50% 
of the luminosity (see e.g lChurch & B alucinsk a-Churchl 1995; 
IChurch et al.ll2002T) . Hence, we fitted the spectrum with a black- 
body + power law model. The best fit gave kT ~ 1 keV, and a 
power law slope 1.6± 0.6;^ 2 /dof =7/14; however, the blackbody 
component was not well constrained, and only contributed ~5% 
of the flux, significantly less than expected for a neutron star sys- 
tem. We infer from the lack of a strong thermal component that 
the primary is more likely to be a black hole than a neutron star; 
however, we cannot rule out a neutron star primary. 

3.2. The X-ray lightcurve 

We obtained the conversions from intensity to 0.5-10 keV flux 
for each Swift observation and the 2010 June Chandra observa- 
tion, assuming a power law with photon index 1.6, with A^h = 
4.2xl0 21 atom crrT 2 . The conversion factor varied by up to 25% 
between observations, due to changes in off-axis angle. 

We present the lightcurve of XMMU J004215.8+41 1924 
covering the outburst from 2010 May 27 to 2010 July 20 in 
Fig. |2] Each point is plotted at the midpoint of the observation; 
the points from the June 18 and June 24 Swift observations, and 
the July 20 Chandra observation, are 3cr upper limits. The 2010 
July Chandra observation yielded a 3<x upper limit of 2xl0 36 
erg s~\ suggesting that XMMU J004215. 8+41 1924 may have 
returned to its pre-outburst level. 

The deepest observation made of XMMU 
J004215. 8+41 1924 is Chandra Obs 1575, with a -40 ks 
duration. The 3<x upper luminosity limit is 6xl0 35 erg s _1 . 
Hence, the variation in luminosity is a factor ^200 between the 
outburst peak and quiescence. 



XMMU J004215.8+41 1924 was observed by Chandra three 
times in 2010 prior to the detection of the outburst, in January, 
February and March; no observation was made in April, as M3 1 
was behind the sun. We found no strong detections in any of 
these observations; the 3cr upper limit for the March observation 
was ~2xl0 36 erg s _1 . Hence the 2010 outburst lasted at least 30 
days, but not more than 140 days. The outburst appears to be of 
similar duration to the one in 2006. 

We have found no evidence of additional outbursts in the 
120 Chandra observations and 90 Swift observations of XMMU 
J004215.8+411924. 

3.3. The X-ray position of XMMU J00421 5.8+41 1924 

For the 2006 July Chandra observation, we chose 12 bright X- 
ray sources associated with GCs for the registration. Five GCs 
had unacceptably large uncertainties in their X-ray positions 
(>0.2"), and these were removed from the registration process. 
Additionally, one of the GCs showed an unusually large dis- 
crepency between the X-ray and optical positions, so it was dis- 
carded and a new solution was found; the final r.m.s. offset be- 
tween Chandra and LGGS was 0.07" in RA and 0.12" in Dec. 
Our best registration solution yielded RA = 00:42:16.063 Dec 
= +41:19:26.73, with 0.17" uncertainty in RA and 0.2" uncer- 
tainty in Dec. Combining these with the position uncertainties of 
the LGGS images gives an error circle with 0.3" radius. 

Our final registration of the 2010 May Chandra observation 
utilised five GCs; the final r.m.s. offset between Chandra and 
LGGS was 0.19" in RA and 0.16" in Dec. Our best location in 
this observation was RA = 00:42:16.037, Dec = +41:19:26.63, 
with 0.17" uncertainty in RA and 0.28" uncertainty in Dec. 
Combining these with the LGGS uncertainties results in a 0.3" X 
0.4" ellipse. The best fit positions in the two observations are off- 
set by 0.3". Therefore it is most likely that both X-ray outbursts 
come from the same source. 

3.4. The Einstein outburst 

We note that our 2010 position is 1.6" from an X-ray source 
detect ed in the Einstein H RI ~30 years previously (source 
17 in Crampton et al. 1984); the uncertainty in Einstein po- 
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Fig. 3. Details of HST images of the environs of XMMU 
J004215.8+41 1924 from 2006, August, i.e. during outburst (left) 
and 2004, October, in quiescence (right). Uncertainties in the X- 
ray position are represented by a circle for the 2006 Chandra 
observation and by an ellipse for the 2010 May Chandra obser- 
vation. 

sition is -10", and includes no known X-ray sources other 
tha n XMMU J0042 15 .8+41 192 This source was not detected 
by lTrinchieri~& Fabbiano| (Il99ll) . who used a source detec- 
tion radius of 6.7" for the Einstein HRC; however, XMMU 
J004215.8+41 192 was observed at a high off-axis angle, and the 
photons were spread over a wider area than the detection cell. 

Extracting a circle with 20" radius around the position of 
XMMU J004215.8+41192 yields 76 counts, while a nearby 
source free region of the same size yields 47 counts. The net 
exposure time for the observation is 28564 s. We calculated a 
conversion from 0.2-4.0 keV intensity in the HRC to 0.5-10keV 
unabsorbed flux using WebPIMMS: 1 count s -1 = 7.1xl0 10 
erg cirr 2 s _1 . Hence we estimate the 0.5-10 keV luminosity of 
XMMU J004215.8+41 192 to have been ~ 5 x 10 37 erg s" 1 in the 
1979 January 13 observation. 

3.5. The search for counterparts 

The 2006, August HST observation was registered to the LGGS 
using seven stars, with a r.m.s offset of 0.02" in both RA and 
Dec. The 2005, May HST observations was registered with 5 
stars, resulting in r.m.s. offsets of 0.019" and 0.05" in RA and 
Dec respectively. 

In Fig [3] we present details of the HST images from 2006, 
August (left) and 2005, May (right). These images are super- 
posed with a circle that represents the X-ray position from the 
2006, July Chandra observation, and an ellipse to represent the 
2010, May position. We find no compelling evidence for a coun- 
terpart during outburst down to rae ~25.5, or my ~26, let alone 
evidence for variability that would confirm the association. 

At a distance of 780 kpc, M31 has a distance modulus 
of 24.45 magnitudes. The observed column density is signifi- 
cantly higher than the Galactic line-of-sight absorption (6xl0 20 
H atom cm" 2 ); us i ng th e empirical relationship obtained by 
Pred ehl & SchmitJ (1 19951) . we can expect ~2 magnitudes of V 
band extinction, and ~3 magnitudes in B band. Hence we place 
upper limits on the counterpart of Mb ^ -2, and My ^ -0.5. 

4. Discussion 

Two mechanisms could be responsible for the huge variation in 
mass accretion that resulted in the observed outburst (see e.g. 
IWilliams etaD 120061 for a discussion of M31 transients). The 



system could be a HMXB with a long, eccentric orbital period 
where accretion is i ntensified nea r peria stron; this scenario was 
initially favoured by lHaberl et al.l d2006l) after mistakenly iden- 
tifying a counterpart in the Swift UVOT image. Alternatively, 
the system could be a low mass X-ray binary (LMXB) with an 
unstable accretion disc that oscillates between a co ld state (qui- 
escenc e) and a hot, ionised state (outburst), see e.g. Dubus et al. 
d200lL and references within). Here we consider the observa- 
tional constraints on both scenarios. 

4. 1 . Constraints on a HMXB system 

The known counterparts of HMXBs in the SMC have apparent 
V magnitudes in the range 13 ^ my 5 18, and B — V in the range 
-0.32 < B - V < 0.06 (see e.g. ICoe et alJl2005t lAntoniou et all 
2009). For a distance of ~60 kpc, this equates to -6 5= My 5 
— 1, all brighter than our threshold of My ^-0.5. It is therefore 
unlikely that a Be star is hidden by the local absorption. 

Furthermore, variations in accretion rate on the orbital cycle 
are of course periodic, and we see no evidence for other out- 
bursts in our ~ 120 other Chandra monitoring observations. Since 
the outbursts lasted at least ~30 days, and the frequency of our 
monitoring is once per ~30 days, we would expect coverage of 
other outbursts. 

If the outbursts in 1979, 2006 and 2010 were due to peri- 
odic accretion near perihelion, then a whole number of orbital 
cycles every ~1430 days would be required. Hence, XMMU 
J004215.8+41192 would have had an unabsorbed 0.5-10 keV 
luminosity ~10 38 erg s " 1 during the 2006 , July 2 XMM-Newton 
observation. However. fVoss et alj d2008h did not detect XMMU 
J004215.8+41 19 during that observation. We therefore conclude 
that XMMU J004215.8+41 1924 is likely to be a LMXB. 

4.2. Constraints on a LMXB system 

The optical emission of X-ray bright LMXBs is dominated by 
the accretion dis c; since larger discs exist in systems with longer 
orbital periods, Ivan Para diis & McClintock d 19941) derived an 

empirical relation between the optical luminosity (Ly), the X- 

i ii 

ray luminosity (Ly) and the disc radius (R) : Ly cc L x R. The 
range in absolute magnitudes of the sample was -5 5 My 5 5. 

Including the relation between R and orbital period, and 
defining 2 as (L x /L Edd ) 1/2 (f/lhr) 2/3 , they further found that the 
relation 

My = 1.57 + 0.24 -2.27 (+0.32) logS (1) 

produced good results over three orders of magnitude in E. This 
relation allows us to estimate the range in orbital period for 
XMMU J004215.8+411924. Due to the high local absorption, 
the upper limit to the orbital period is not particularly constrain- 
ing: 540 hr for a neutron star and 5 130 hr for a black hole, 
longer than typical LMXB periods. 

The X-ray spectrum from Observation 11838 is best fit- 
ted by a power law fit with spectral index ~1.6; this is seen 
in neutron star and black hole binaries at low accre tion rates 
(Ivan der Klisl [19 94). However. iGladstone et al.l d2007l) showed 
that neutron star LMXBs don't exhibit this behaviour at 0.01- 
1000 keV luminosities ^10% of the Eddington limit. Neutron 
star LMXBs at higher luminosities have two component emis- 
sion spectr a, with a thermal c omponent contributing -10-50% 
of the flux dChurch et al.l2002l) . Since the 0.5-10 keV luminosity 
of XMMU J004215.8+41 1924 is -70% of the Eddington limit 
for a 1 .4 M neutron star, the primary is likely to be a black hole 
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dBarnard et al.l l2Q(58h : this conclusion is supported by the lack 
of a strong thermal component in the two component model. 
However, the quality of the data prevents us from excluding a 
neutron star primary. 
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